Single Spin Superconductivity: 
Formulation and Ginzburg-Landau Theory 



Robert E. RudcU and Warren E. Pickett* 
Naval Research Laboratory, Washington, DC 20375-5345 
t SFA, Inc, 1401 McCormick Drive, Largo MD 20774 
(February 7, 2008) 



Abstract 



We describe a novel superconducting phase that arises due to a pairing in- 
stability of the half-metallic antiferromagnetic (HM AFM) normal state. This 
single spin superconducting (SSS) phase contains broken time reversal sym- 
metry in addition to broken gauge symmetry, the former due to the underlying 
magnetic order in the normal state. A classification of normal state symme- 
tries leads to the conclusion that the HM AFM normal phase whose point 
group contains the inversion operator contains the least symmetry possible 
which still allows for a zero momentum pairing instability. The Ginzburg- 
Landau free energy for the superconducting order parameter is constructed 
consistent with the symmetry of the normal phase, electromagnetic gauge in- 
variance and the crystallographic point group symmetry including inversion. 
For cubic, hexagonal and tetragonal point groups, the possible symmetries of 
the superconducting phase are classified, and the free energy is used to con- 
struct a generalized phase diagram. We identify the leading candidate out of 
the possible SSS phases for each point group. The symmetry of the super- 
conducting phase is used to determine the cases where the gap function has 
generic zeros (point or line nodes) on the Fermi surface. Such nodes always 
occur, hence thermodynamic properties will have power-law behavior at low 
temperature. 
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I. INTRODUCTION 



The pairing theory IJTJJ^] of superconductivity and superfluidity is based on a normal 
state with time reversal symmetry and inversion symmetry. The former symmetry requires 
that the two spin directions are related by symmetry, and spin-space rotations have played a 
central role in the classification of broken symmetry phases. The latter symmetry is sufficient 
to ensure that states with identical spin directions at k and —k are degenerate. Electrons 
at k and —k can then be paired to have zero total momentum and can be classified by spin 
state as singlet or triplet. The question of inversion symmetry and time- reversal-breaking is 
central in the specification of the pairing states, |0HTD[ and in this paper we present some 
new aspects of these relationships. 



It has recently been pointed out |TT[ that there is a normal state with broken time 
reversal symmetry that has a pairing instability in direct analogy with that in BCS theory. 
||T|| This magnetically ordered normal state, in which inequivalent up-spin and down-spin 
magnetizations cancel exactly, has been termed "half-metallic antiferromagnet" (HM AFM) 
by van Leuken and de Groot. Jl2| This normal state, which we describe in Sec. [II], has 
considerable theoretical and probable technological interest in itself. Its lack of macroscopic 
magnetization means that considerations of pairing do not have to confront the question of 
competition between superconducting order and a pre-existing magnetic field. 

In this paper we outline in more detail the characteristics of this "single spin supercon- 
ductivity" (SSS) phase. In Sec. |TJ we review the characteristics of the HM AFM state, which 
is the precursor normal state of the SSS. In Sec. [HJ we show that the phenomenon maps 
onto the BCS model with minor but non-trivial changes. Comparison to liquid 3 He, con- 
ventional BCS superconductors, and more exotic heavy fermion superconductors in Sec. ^ 
demonstrates that the HM AFM state has the minimum symmetry required of the normal 
state to allow q = pairing instability, and that in SSS theory inversion symmetry plays 
a role analogous to that of time-reversal symmetry in BCS theory. In Sec. [V| we provide 
the classification of all possible order parameter symmetries for high symmetry crystal point 
groups, and present results of a symmetry analysis of the Ginzburg-Landau free energy that 
enumerates the possible SSS states for cubic, tetragonal, and hexagonal crystals. 



II. THE NORMAL STATE 

A. Half- metallic ferromagnetism 

A half-metallic (HM) ferromagnetic (FM) electronic structure arises in a ferromagnetic 
material when the Fermi level (Ep) of one spin direction lies within a gap in the spectrum 
of the other spin direction. |l3|Jl^| The gap may occur in either the majority or the minority 
channel. In either case, we will take the up channel to be the metallic one. We also confine 
this discussion to stoichiometric compounds, which have an integer number of electrons per 
cell. The system of up spins then forms a metallic fermion liquid, while the down spins form 
an insulating system that may be thought of as an inert background for the purposes of 
studying low temperature, low energy processes. This specific occurrence, the placement of 
the Fermi level of the metallic channel in a gap of the other channel, defines half-metallic 
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character: the up-spin "half" of the electrons is metallic, while the down-spin "half" is 
insulating. Fig. |l|(a) shows a model spectrum of exchange split bands that leads to HM 
character. 

Half-metallicity leads to several features of a crystalline solid that are qualitatively dis- 
tinct from conventional metallic ferromagnets. Unlike in a conventional FM, electron trans- 
port is 100% polarized, and there are no allowed low energy spin flips. In a common FM, the 
spin moment is a continuous quantity whose value is determined by the balance of exchange 
energy and kinetic energy. In a HM FM, however, the spin moment is constrained to be 
precisely an integer A4. This is so because the insulating channel contains an integer number 
of filled bands, and hence an integer spin A/j, the cell contains an integer number of elec- 
trons J\ftot, so the metallic up channel contains an integer number of spins A/j = Aftot — A/j.. 
Then M. = A/f — A/j. is an integer, and the moment is A4//B per cell. This follows for any 
placement of Ep within the gap in the down-spin density of states. The application of a 
magnetic field H shifts up and downs spin bands by ig/x^H but does not change the occu- 
pation or the net spin moment. Hence the spin susceptibility is precisely zero, which is a 
direct consequence of the lack of low-energy spin flips. There is then no Stoner continuum 
to damp spins waves by single spin flips of carriers. In fact, the situation can be categorized 
as extreme spin-charge separation in the carrier system, in which the spin degree of freedom 
has been separated from the charge fluctuations and frozen out entirely. In this paper we do 
not address possible effects due to spin waves. 



Perhaps the simplest example of a HM FM is CrC>2, in which the moment is 2//b. [|15j de 
Groot and collaborators Jl3| have identified calculationally various Heusler alloys that are 
likely HM FMs, and experimental work on several members (especially UNiSn and NiMnSb) 



has been reported. [16,17] Pickett and Singh [ffq] presented theoretical evidence that the 



colossal magnetoresistance manganites, viz. La2/3Cai/3Mn03, are HM in their low temper- 
ature FM phase. Recently several candidates for HM behavior have been found with the 
double perovskite crystal structure. |0| These examples indicate that half-metallic character 
is not a rare phenomenon. 



B. Half-metallic antiferromagnetism 

It may occur that the integer spin moment M. in a half-metallic system is zero. This 
situation has been termed half-metallic antiferromagnetism (HM AFM). |12] Its properties 



are like that of the HM FM discussed above, with one essential difference: there is no 
macroscopic magnetic field. The HM AFM has 100% polarized charge transport without 
any net magnetization. It must be kept in mind that the HM AFM is not antiferromagnetic 
in the usual sense of the term, as there is no symmetry operation that connects spin up 
and spin down states or densities. In fact, it is essential that the two spins channels are 
electronically (and thus chemically) distinct, so a gap can occur in one channel only at the 
same band filling. A model illustrating this situation is shown in Fig. |l](b). This model 
uses the same two-band form of Fig. |l](a), but the bands for the two spin channels must be 
displaced in opposite directions, reflecting the necessary inequivalence of the channels. 

In a HM AFM the spins are precisely balanced, so there is no majority or minority spin. 
In this paper we will call (when the need arises) the metallic channel the "up" spin and 
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the insulating channel the "down" spin. For low energy and low temperature processes the 
insulating channel becomes inert and drops out of consideration, van Leuken and de Groot 
|T2|| have suggested one quintinary compound, intermetallic VyMnFegSbyln in a Heusler-like 



crystal structure, that should be a HM AFM. One of the present authors [|ll] , |19| has found 
candidates for HM AFM states within the class of magnetic double perovskites. An example 
of a HM AFM spectrum is shown in Fig. 1 for the double perovskite compound La2VCu06, 
calculated using accurate spin density functional methods. [JT3] This is nominally a Cu 2+ , 
V 4+ compound for which both ions have spin |. When the spins are parallel (top panel) the 
spectrum is that of a conventional metallic FM. When the spins are antialigned, however, 
the Fermi level falls where the spin up density of states is large, but within a gap in the spin 
down channel. This compound, and other double perovskites, are discussed in more detail 
elsewhere ||19|| . 



We leave further discussion of proposed HM materials to future papers, and address 
below the pairing instability of the HM AFM normal state, and its consequences. 



C. Experimental Consequences 

An important practical consideration is how a HM magnetic material can be identi- 
fied. The anticipated properties (questions of many-body corrections [^] aside) have not 
previously been enumerated. We include a partial list here to provide guidelines. The gen- 
eral feature of course that as the temperature is lowered through the magnetic ordering 
(Curie or Neel) temperature Tm, the material changes from an nonmagnetic (conducting or 
non-conducting) system to a metallic magnetically ordered system at low T where the spin 
excitations are frozen out. 

a. Metal with fully polarized transport at low T. Metallic resistivity, but vanishing magne- 
toresistance at low T. There is no clear signature in the Hall or Seebeck coefficients. At 
intermediate temperature there may be a negative magnetoresistance for a HM FM, 
reflecting the field induced increase in magnetic order and reduced spin scattering as 
the carriers in one channel become non-conducting. For HM AFM, this field induced 
effect will not apply. 

b. Magnetic order. There is no obvious signature of the HM character in the spin wave 
spectrum or temperature dependence of the magnetic order parameter. 

c. Vanishing spin susceptibility. However, core diamagnetism, van Vleck (orbital) param- 
agnetism of the metallic channel, Landau diamagnetism of the insulating channel, and 
temperature variation of the net order of the local moments will make the magnetic 
susceptibility difficult to analyze. 

d. Non-Korringa behavior in NMR. This technique may provide the most direct indication 
of HM character. The longitudinal relaxation rate Tj -1 , which is a measure of the 
conduction electron spin flips, is proportional to the product of the densities of states 
of each spin channel [Nj(E_f)Nl(E,f)], which vanishes for a HM phase. The Knight 
shift, normally dominated by the spin susceptibility in normal metals, should be small. 
An NMR study of the proposed HM magnet UNiSn has been reported. p2i 
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Spin-polarized electron spectroscopies. At first glance, these spectroscopies seem ideal, 
but both photoelectron emission and STM tunneling are sensitive to surface properties. 
In addition, the magnetic order may be different at the surface, mitigating against HM 
character in the surface region. Spin-polarized photoemission studies of Cr0 2 were 
inconclusive. [ f23| . |15" [ 



Spin-polarized positron annihilation. This technique, which takes advantage of the 
natural polarization of the positron beam, has been claimed to establish within narrow 



bounds that NiMnSb is a HM ferromagnet. pi It is a bulk probe. 



g. Thermodynamic properties. We will show that SSSs necessarily have point or line 
nodes of the gap (not always the case in previously studied cases of triplet pairing). 
The resulting gapless excitations lead to heat capacity, penetration depth, thermal 
conductivity, etc. that have power-law in T (or uj) rather than exponential. 

h. Tunnelling. Tunnelling between an SSS and a ferromagnet will show a strong de- 
pendence on the direction of magnetization of the ferromagnet. Josephson coupling 
between an SSS and a singlet pairing superconductor should not occur. 



III. PAIRING INSTABILITY IN THE HM AFM 

The HM AFM is a single component fermion liquid as a result of underlying magnetic 
order and electronic structure that renders one spin channel insulating. The lack of a micro- 
scopic magnetic field in a HM AFM suggests the possibility that a superconducting instability 
may occur in the metallic channel. The Cooper instability [0 is spin-blind: the two fermions 
that undergo the pairing instability can have antiparallel spins as in BCS theory, or they can 
have parallel spins as in 3 He, and the instability is straightforwardly extended to a spinless 
fermion system. We now show that this instability maps directly onto the BCS model of 
superconductivity ]1[ in a simple but not quite trivial manner. This new superconducting 
state has been called single spin superconductivity fLlf . 



A. Formal Relationship to BCS Theory 

In BCS theory an electron in state K is paired with its time-reversed partner TK. 
K = (k, |) is an index that together with its partner TK = (—k, j) exhausts all states on 
the (| and I) Fermi surface(s). In a SSS, an electron in state K = (k, +) is paired with its 
inversion partner TK = (—k, — ). To cover all states on the Fermi surface once only, K must 
range over only half of the Fermi surface, say the 'top' half with k z > (hence the notation 
'+'), and states with k z = can be assigned to '+' and '— ' components of the pairs also. 

In terms of the general two-body interaction 

all terms except ki = —k 2 = k, k 3 = —k^ = k' are irrelevant for pairing. To count the pair 
states properly, the full Brillouin sums in the above expression must be expressed in terms 
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of pair indices K. Using the anticommutation relations of the electron field operators at, 
and defining the pair annihilation operator = a^a_^, the resulting interaction is (making 
the notational simplification V^_g^, _p — > V^g,), 

v P air = E E [%, + ^ - - v_ u ] b% = EE ^ A- (3-2) 

fe fe' fc fe' 

The '+' sign indicates that the summation extends only over the 'top' half k z > of the 
Brillouin zone. The symmetry of the matrix element noted by Sigrist and Ueda |2j results 
in all terms in the braces being identical, so U^p = 4Vgg ; . We return to the implications 

of the form of the matrix element in the next subsection. From V^p = (k 1 , —k'\V\k, —k) = 
(K' ,XK'\V\K,XK) = Vk,k>, and similarly for Uk,k', and the expression for the kinetic en- 
ergy in terms of pair labels (running over half of the Fermi surface), the pairing Hamiltonian 
then is 

Hp a ir = E e k( a K a K + aj K a IK ) + E E U K,K'b ] K b K i (3.3) 

K K K' 

where bx = cikclik- Note that we do not use the BCS convention of separating out a negative 
sign from the interaction matrix elements. 

The Bogoliubov-Valatin transformation [f24j is analogous to its form in BCS theory, 

a K = u K a K + v K a\ K , (3.4) 

a iK = u iK a iK + Vj K a K , (3.5) 

and the formalism of SSS pairing theory maps onto BCS theory. Specifically, the SSS ground 
state is 

% = X\{u K + v K b ] K )§ vac (3.6) 

K 

and the gap function is given by 

A K = J2U K ,K'(b K }. (3.7) 

K' 

The gap is a scalar, i.e. it has no spinor indices, which distinguishes SSS from all supercon- 
ducting systems discussed previously. The gap equation is formally identical to BCS: 

A « = ~ E SF A ^' tanh (¥ E k) (3-8) 

K , Z£jk< v 7 

where (3 = 1/kT is the inverse temperature. With \i denoting the chemical potential, 

E K = Ei=\{ H -tf + \±tf\\ (3.9) 

is the quasiparticle excitation energy, which is even in k. 

It is useful to express the gap equation as usual in terms of a full Brillouin zone summa- 
tion. It is readily verified that extending the sum over the full zone introduces an expected 
factor of |, and the gap equation, 

A fc = ~E T^Ag, tanh (IfiEp) (3.10) 

where Wr.p = \U^p = ^kk" indicates the final formal equivalence to the BCS equation. 
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B. Simple Consequences of Single Spin Pairing 



The combination of matrix elements in Eq. (|3.2| ), and the symmetry noted by Sigrist and 
Ueda, @ 

= -V-ifi = ~ V k,-k> = V-t-fr (3- 11 ) 

indicate explicitly the /c-space structure that is necessary for SSS pairing. A fc-independent 
attractive potential V^p = —V, which leads to singlet pairing in the BCS model, contributes 

nothing toward SSS pairing; likewise, a /c-independent repulsion is harmless. The simplest 
form of such coupling, which is odd in both k and k', is of the form Wk ■ k . When this is 
substituted into the gap equation Eq. ( |3.10| ) (see below), a non- vanishing solution requires 



W < 0. Thus the pairing interaction must be attractive for small angle ("forward") pair 
scattering and repulsive for large angle scattering of pairs. This behavior is reminiscent 
of the situation in high T c theory, where the spin-fluctuation picture has an everywhere- 
positive interaction, which peaks at large q (more specifically, at (it /a, rr/a)). |25| That type 
of interaction favors a d x 2_ y 2 symmetry of A for singlet pairing. 
Using the pairing interaction in normalized form 



W W = -\W\^ (3 ' 12) 



where kp is the Fermi wave vector, the T=0 gap equation can be solved readily for several 
trial gap symmetries. Gap functions of the form Ag ex d ■ k, for some constant vector d, 
give the lowest order possibilities (in terms of polynomials of the components of k). If d is 
real, or purely imaginary, its direction can be taken as the z axis (d= (0,0,1)) so Ag oc k z 
with a line of nodes on the equator of the Fermi surface. Complex d can be represented by 
d=(l,i,0)/y/2, in which case Ag cx k x + ik y with point nodes at the poles. We also consider 
the "highest symmetry" single dimensional function t[ oc k x k y k x {k^ — k y ^{ky — k 2 z ){k 2 z — k'l) 
(see the following Section) as an exotic possibility - it has nine lines of nodes. 
We treat the usual weak-coupling case, where 

E ^ « ^(0) «, (3,3) 



4-7T J-w c J Air 

k 

where N(e) is the density of states per spin which is assumed to be constant over the energy 
scale ui c of the pairing interaction. We display in Fig. 3 the resulting T=0 gap value, relative 
to the energy cutoff u c , versus the coupling strength A = iV(0)|W| for these gap symmetries. 
The BCS result is given for comparison. Given the same coupling A, it is evident that the 
zero temperature gap magnitude is comparable to the BCS value, even for the rj~^ function. 

It is straightforward to obtain the limiting behaviors of the zero temperature gap from 
Eq. ( |3.10D . In the weak-coupling limit, 



A rms = 2^ c e- 1 / A C 1 (l + C 2 e" 2 / A ), (3.14) 
while for large A the asymptotic form is 
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D 2 

A rms = u c X(D 1 - — ). (3.15) 

Here C±, C 2 , D x , D 2 are symmetry dependent constants (c/. Appendix B). This latter relation 
explains the linear behavior at larger A that is evident in Fig. 3. Note that strong coupling 
corrections will change this to a v^A behavior. |26 



At T = T c , Et — > eg (we take /x = 0). Taking k along d (or along one non-zero component 
if d is complex), the only difference from the corresponding BCS equation is an angular inte- 
gral. This integral is unity, however, due to the expansion of Wt p in normalized functions. 
Thus the equation for T c versus coupling constant A = iV(0)|W| is identical in form to that 
of BCS: 

UfiM, (3.16) 
A Jo e 

where (3 C = l/fcsT c . 

Although there is every reason to expect that SSS will arise in the appropriate transition 
metal (or / electron) compounds, there is not yet any expectation of high T c . For one thing, 
the necessary interaction Wk,k' is of a particular kind (see above); however, this is also the 
case for the d wave scenario in high T c cuprates. More to the point, however, is that at 
larger T transverse spin fluctuations increase strongly and tend to reduce the AFM order 
parameter, finally causing the system to revert to the paramagnetic state above the Neel 
temperature. The manner in which the HM AFM state will thereby be weakened has not 
yet been explored. 



IV. SYMMETRY: RELATIONS TO PREVIOUS THEORY 

It is instructive to clarify the relationship between the degree of symmetry of the normal 
state and the degree of richness of broken symmetry in the condensed phase. 3 He has the 
highest symmetry possible in its (liquid) normal state. It has continuous real space C and spin 
rotation S symmetries, it has time reversal T and inversion X, and of course gauge symmetry 
U(l). The group of its normal phase then iszCxiSxXxXxU(l). In the condensed superfluid 
phase, U(l) is a broken symmetry and one or more of the other symmetries can be broken 
concomitantly. Much work has been done to characterize the more likely cases among the 
infinite possibilities (infinite because the relative angular momentum L of the pair can be 
any non- negative integer). The observed phases correspond to particular states within the 
(pair spin and orbital angular momentum quantum numbers) S = 1, L = 1 complex that is 
described by an 18 component order parameter in general. 

When considering pairing in a crystal, the continuous real space rotation symmetry is 
replaced by the finite crystalline point group G (see Table Q, where the classifications of this 
subsection are collected). The group of the normal state then is GxS x T x XxU(l). The 
number of broken states is finite because the space of basis functions of irreps is spanned by 
a few small-L sets. The necessary values of L for cubic crystals are presented under the BCS 
case in Table |. The allowed symmetries in cubic, tetragonal, and hexagonal crystals have 
been exhaustively categorized. f2|]27|| 
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In crystals with strong spin-orbit coupling, it may be appropriate to consider the spin 
as frozen into the crystalline lattice, in which case spin rotation is no longer a separate 
symmetry. (On this matter there are arguments both pro and con in the literature.) Then 
the normal state symmetry is lowered further, and the number of distinct broken symmetry 
states is further reduced. The next lowest symmetry situation, where inversion symmetry is 
absent from the point group, has been considered by Poluektov [28] with the spins frozen to 
the lattice, distinction between singlet and triplet pairing vanishes and the order parameter 
becomes a low symmetry combination. 

For the SSS discussed in this paper, the two spin systems are inequivalent, so time-reversal 
and spin rotation symmetries are strongly broken by the normal state. The symmetry 
of the normal state is described by GxXxU(l). The number of distinct possibilities for 
broken symmetry states, which we enumerate below, is reduced still further. In fact, this 
state of affairs gives the lowest possible symmetry normal state that still allows pairing in 
the usual sense of zero net momentum Q of the pair. Inversion symmetry ensures that 
e£ = e K = txK = so that if k lies on the Fermi surface then so does — k, and these two 

states can pair to total momentum Q = k + (—k) = 0. Without inversion symmetry, z.e.for 
GxU(l) normal state symmetry, Q = pairing is not allowed. This result indicates that 
inversion, not time-reversal symmetry, is the minimal symmetry requirement for a pairing 
instability. For convenience, inversion I will be considered part of the point group G below. 

Because for a SSS state time-reversal symmetry is already broken in the normal state, 
it is natural to expect (and indeed we find) that broken symmetry phases with unusual 
properties are likely to arise. Such considerations occupy the rest of the work reported here. 



V. GINZBURG-LANDAU FREE ENERGY 
A. Allowed Superconducting Phase Symmetries 

The spontaneous symmetry breaking at the superconducting phase transition is governed 
by the free energy. Below T c the minimum is a superconducting phase whose order parameter 
breaks the U(l) gauge invariance and possibly other symmetries as well. This physics may 
be described through a phenomenological Ginzburg-Landau free energy that describes the 
mean field theory of the superconductor in terms of a few parameters related to matrix 
elements of the effective potential. Despite its simplicity, the mean field free energy captures 
all of the generic information about the allowed symmetries of the superconducting phase. 

The order parameter describing the Cooper pair condensate is taken to be the gap func- 
tion, 

^5 = J2W lP ( H ,a_ p ), (5.1) 

k' 

which was introduced in Equation ( |3.7p . Below T c the gap function is non-zero, and it 
transforms under the full symmetry group of the normal phase. In particular, it is not 
invariant under the U(l) Abelian gauge symmetry (and the concomitant global U(l) of 
electron number), since Cooper pairs have electric charge —2 (electron number 2). The U(l) 



9 



symmetry is broken to the cyclic group Z2. The gap function may break other symmetries 
as well, and this is the main focus of this Section. 

The gap function must form an irreducible representation (irrep) of the symmetry group, 
Q = G x U(l), at least in the vicinity of the phase transition. This follows from the fact 
that near the transition, the gap function satisfies the linearized gap equation, which for a 
spherical Fermi surface is given by 

A k = - 7 (o;e/*r c ) Jdk'W^Ay 

tanh ( 5 - 2 ) 



with 7(0;) = N(0) / dy 



y 



where ui c is the cutoff for the interaction. This is a (/-invariant eigenvalue equation where A 
is an eigenvector and therefore must transform as a member of an irrep of Q. As usual, we 
start with the assumption that the transition temperature for the first irrep to condense is 
much higher than that of the others. This lets us focus on each irrep separately. 

Under a U(l) gauge transformation —>■ e~ vp at 1 A — ► e~ 2iv A. A may transform under 
the point group, G, as well. Consider the case where there is an element of G that sends 
A to another function which is not related to A by a gauge transformation; i.e. there is 
g G G such that gA = A' where |A| 7^ |A'|. This breaks the point group symmetry to a 
subgroup H of G that does leave A invariant up to a gauge transformation. This non-trivial 
form of spontaneous symmetry breaking occurs exactly when A is in an irrep r of G whose 
dimension is greater than one |29 . 



The overall symmetry breaking scheme may be described as follows. When A is in 
the irrep T of G, there is a maximal subgroup H of G under which A transforms as a 
one dimensional irrep T'. Then Q is broken to the little group H, Q — > H. All such 
decompositions of irreps of the crystallographic points groups are tabulated in "compatibility 



tables" ||30|| . It is merely a matter of looking up the maximal subgroup H in which each 
distinct one dimensional irrep appears. This classifies all of the possible superconducting 
phases according to symmetry. 

There is one subtlety in this analysis. The gap function is not gauge invariant, so it is 
not a physical observable. The physical residual symmetry group H p h ys may be larger than 
H. For example, the spectrum of quasiparticle excitations Q3.9| ) has the symmetry of the 
gauge invariant product A*A. If A is in a complex irrep, there may exist an element g of Q 
that switches V and its complex conjugate V ; that is, gip = ip for ip G V. Then g leaves |A| 2 
invariant, but it is not an element of H. The residual symmetry of the physical observables 
below T c is the group H p h ys generated by H and g (H p h ys — H x Z 2 ). On the other hand, if 
there is no such g such as when the irrep is real, then H p h ys = H . 

It is convenient to expand the gap function in terms of explicit representatives of the 
irrep V. A general form of the mode expansion of A^ is 

^^EEtfi™')^'™;^^) ( 5 - 3 ) 

r,m fit 

where r] m (T; rij) is a complex coefficient, m — 1, • - • , dimT labels the components of V and 
the indices distinguish polynomials of different degrees. (For multiple or aspherical Fermi 
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surfaces Allen's Fermi surface harmonics |3T| would be used.) A suitable choice of basis 
elements ip(T, m; nj = 0; k) is given in Table PJ A complete basis of higher modes of the 
linearized gap equation may be constructed using standard techniques. 



B. Construction of the Free Energy 

At this point we have identified the possible residual symmetries of the superconducting 
phase. It remains to show how each is realized as a minimum of the free energy. The 
Ginzburg-Landau free energy, F, is a functional of the mean field A^ in which the fluctuations 



about the mean field have been integrated out It must be invariant under the full 

symmetry group Q which describes the physics of the normal state. This constrains the 
combinations of the modes of A^ (|5.3|) that enter F. Only invariant (Tf) combinations 
contribute. In a perturbative expansion of F in terms of Ar, the form of the low order terms 
is highly constrained by (/-invariance. 

The group symmetry imposes a very restricted form. Many terms that could appear 
in the free energy vanish. This is usually computed with a Clebsch-Gordon decomposition 
HJ30|], but we develop a more powerful technique in which the group symmetry is used 
directly. The result is an explicit computation of the allowed terms in the free energy to 
arbitrarily high order in perturbation theory. 

The group symmetry imposes a number of constraints. Gauge invariance requires that 
the polynomial have equal numbers of 77's and 77's. The restrictions due to the point group 
are implemented as follows. Each group operation may be considered to act on the 77's in a 
way that leaves the free energy invariant. If A is transformed under an operation from Q, 
it may be restored to its original form by a linear transformation of the constants rj m ; that 
is, they form the contragredient representation of T of G with charge 2 under global U(l). 
The free energy is invariant under Q, so it must be an invariant polynomial in rj m under 
the action of Q . In addition to the symmetry constraints, the free energy must be real and 
bounded below for stability. 

Invariant polynomials have been studied extensively in the mathematics literature. In 
particular, invariant polynomials for the symmetric and alternating groups |34| and for 



Abelian and non-Abelian gauge groups have been constructed explicitly. These polyno- 



mials play an important role in gauge theory |36| , |37| . Invariant polynomials for G x £7(1) 



require an extension of this theory, and since it has not been discussed in the condensed 
matter literature, we will give some details. 

Since G is a finite group, its action on the coefficients rjj is isomorphic to a direct product 
of finite simple groups-in particular the symmetric groups, S3 and 5*2, the alternating groups 
A3 and A2 and the cyclic groups, Z3 and Z2. The form of invariant polynomials for each of 
these groups is well known, and our task is to form combinations of them that are invariant 
under G x U(l). 

The free energy is almost trivially constructed for the one dimensional representations. 
Gauge invariance requires that the free energy be a function of \i]i\ 2 . This is also invariant 
under the point group operations. The perturbative expansion of the free energy takes the 
form 

F = a|?7i| 2 + /?|?7i| 4 + 7|r7i| 6 H for dimT = 1 (5.4) 
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where a, (3 and 7 are parameters describing the expansion of the effective potential for 
the order parameter in terms of the basis elements of the irreducible representation. Note 
that the k integrals for matrix elements of the effective potential Wzti have been included 
in these parameters, so they encode the physics. In general, there are exponentially small 
corrections to this perturbation series of the form P(|r/i| 2 ) e _1 / ( - Q '' ,?1 ' - 1 where P is a polynomial. 
These non-perturbative corrections can be important at very low temperatures or in strongly 
coupled systems, but they are beyond the scope of this paper. 

The treatment of higher dimensional representations is more involved. The first step 
is to find a basis for the point group irrep that respects the simple group decomposition. 
Each simple group generator should correspond to a specific element of G. This makes the 
symmetries of the invariant polynomials manifest, permitting a straightforward description. 
Using invariant polynomial techniques, we will construct the free energy for and rj 5 of 
O h , 6 of D 6h and Tg of D ih . 

Consider the two dimensional representation of Oh- The standard real basis shown in 
Table [II] does not respect the simple group decomposition. The appropriate basis is composed 
of the complex functions 



= [^i(rj) + z^ 2 (r 3 )]/ v / 2 = k x k y k z (k 
ip' 2 = [^i{T^) -i^ 2 {Ts)]/V2 = k x k y k z (k 



+ uk 2 x + u 2 k 2 y ) 
* + u 2 kl + uk 2 y ) 



(5.5) 



The action of Oh on this basis is isomorphic to Z3 x Z2 x X, where Z3 = {E,Cs a ,C^}, 
Z 2 = {E, C' 2a \ and X is inversion (another Z 2 ). In addition to U(l), we only need to consider 
the following two point group operations: 



G 



3o 







2a 




(5.6) 



Consider a monomial (rj[) Nl {rf 2 ) N2 (^fi)^ 1 {v^) 1 ^ 2 * n ^ ne perturbative expansion of the free 
energy. Invariance under Oh x U(l) imposes the following constraints 

a. U(l): Ni + N 2 ^N x + N 2 . 

b. C 3a . Nx - N 2 - Fi + N 2 = (mod 3). 

c. C' 2a : invariance under r][ — > —r]' 2 ,r] 2 — > —rf v 

This leads to a basis of invariant polynomials of the form 



1 1 2/3 



and 



(5.7) 



They can be reexpressed in terms of the polynomial generators for T 3 of Oh 



Pi 
Pi 
Ps 
Pi 



4|«| 2 



8TZe 
8 (Wi 



1— 1 ^3 

/ \2 



1 |2 



Xm 



(v'lViY 



(5i 
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where the complete set of invariant basis elements are powers of these four elements, in one 
of the two forms 



p+ _ T>n pp p 

r (m,n,p) ~ r l r 2 r X 

p— pn pP pm p 

r {m,n,p) ~ r \ r 2 r 3 r ^ 



(m,n,p) 1 2 3 /g g\ 



(5.11) 



These invariant polynomials have not been constructed previously. 

Using the invariant basis the most general form of the perturbative expansion of the free 
energy may be written 

F{Oh(T 3 )) = ^2 -^(m,n,p) -^(m,n,p)(^l' %) + ^(m,n,p) ^(m,ri,p) (^1 > 

m,n,p (5.10) 

= a (K| 2 + |^| 2 ) + # (K| 2 + |^| 2 ) 2 + 4[3> \ V ' lV ' 2 \ 2 + ■■■ 

where a = F^ 10 ^, (3[ = F^ 20 y j3' 2 = F^ Q1 y etc. For comparison, Sigrist and Ueda use 
the coefficients f3\ — [3[ + p 2 and /? 2 — fi'i- The minima of this free energy will be studied in 
the next Subsection. 

The next case is the two dimensional representations and Tg of D 6 h- The analysis is 
essentially identical to the case of T 3 of Oh- Again the action of D Q h on these two irreps is 
isomorphic to Z 3 x Z 2 x T, where Z 3 = {E, Cs z , C^}, and Z 2 = {E, C 2a; }. The difference in 
the two cases comes from the way Z 2 xl is embedded in the group; e.g. C2 y (T^ ) = IC2 V (Tq ). 
A complex basis is necessary to make the simple group decomposition manifest 

= [Vi + ^/v^ 

= (k x + iky) I \/2 for 
= k y k z (k x + ik y )(k 2 y -3k 2 x )/V2 for Tg 

= (k x — iky) / \pl for 

= k y k z (k x - ik y )(ky - 3k 2 x )j\pl for T 6 . 

The generators for the invariant polynomials of T§ 6 of are the same as those for Tg of 
Oh (|5.8|), and the free energy may be expressed 

F(D<5h(X 5 6 )) = 22 F( m)njP ) ^ > ( m ,n,p)( 7 ?i> 7 72) + ^(m,n,p) ^(m,ra,p) (^l > %) 

rn,n,p (5.12) 

= a (K| 2 + |^| 2 ) + # (l^'J 2 + l^l 2 ) 2 + 4/3 2 |^ 2 | 2 + ■ • • 

This is identical in form to F(Oh(J?z ))■> so the phase transitions take place at the same 
values of a, Pi, 7$, etc. Of course, the symmetries of the phases are different in the two cases. 
Also, the dependence of the parameters a, A, . . . on physical quantities such as couplings, 
masses, the temperature and the pressure are different, so systems with different normal 
state symmetries do not sit at analogous locations in the superconducting phase diagram in 
general. 

Next we consider the two dimensional representation Tg of The action of D±h on 
this irrep is isomorphic to S2 x Z 2 ; that is, it permutes 771 and 772 and changes their signs. 
In particular, S 2 = {E, C 2e }, and Z 2 = {E, C 2y }, {E, C 2z }. The generators for the invariant 
polynomials of of D Ah are 
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Pi 

Pi 

P 3 
Pi 



1 2 i I |2 

2 i 2 2 /| |2 i I |2 N 2 

V1+V2 - (M + lr?ss 
I^i^l 2 

4 (l^il 2 - M 2 ) 1m [(r/i^ 2 ) 2 



4 [Im (r7ir/ 2 )] 2 



(5.13) 



The basis elements are generated exactly as in the case of T 3 of Oh (|5.9|) , and the free energy 
is given by 

F(D 4h (T^ )) = ^ ^(m,n,p) P^n,n,p)(Vl,V2) + F( m , n ,p) P (m,n,p)(Vl, V2) 



m,n,p 



= a(\ Vl \' 2 + \r] 2 \ 2 ) + (3x (Ir/xl 2 + |^ 2 | 2 ) 2 - 4f3 2 [Im(^ 2 )] 2 + 4/3 3 |^ 2 | 2 + • • • 

(5.14) 

The coefficients agree with those used by Sigrist and Ueda @ to fourth order, except for (3 3 
which differs by a factor of 4 with theirs larger. 

The final irreps are the three dimensional representations and of Oh- The action 
of Oh on each of these irreps is isomorphic to S3 x Z 2 ; that is, it permutes rji, r/ 2 and r] 3 and 
changes any of the signs. Of course, the correspondence between specific group elements and 
these transformations differs in the two cases. In Tj, S3 is generated by O35 and CiCy and 
the three copies of Z 2 are generated by the reflections a x ,a y and a z . On the other hand, 
in , 5*3 is generated by C 3 s and C' 2b and the Z 2 actions are generated by the reflections 
<j y ,<j z and a x . The details of the construction of these invariant polynomials and those for 
the two dimensional irreps are presented elsewhere [Q. The generators for the invariant 
polynomials of Tj 5 of Oh are 



Pi 




m\ 2 


+ M 2 + 


V3\ 2 










(50(6)) 


Pi 






il 2 + \V2V3\ 


2 + 


l™| 2 








(f/(l) 3 x O h ) 


P3 




ViViVs] 2 












(U(l f x O ft ) 


Pi 




2 1 2 1 2 2 

V1+V2+ V3 










(C/(l) x 50(3)) 


A 




2—2 , 2—2 1 
V1V2 + V2V3 + 


2— 


2 2 p2 
1 ~ *2 








P(, 


= Tie 


1 \2 / 2—2 

{V2V3 


- \V2V3 | 2 ) 


+ eye. 






P7 


= Im 


vlvl (M 1 


| 


m\ 2 )' 


■f eye. 






Ps 


= Tie 


1 |2 2—2 / 

\Vi\ V2V3 { 


2^1 


r - 


|2 




+ eye. 


Pg 


= Im 


4—2—2 , 4—2—2 , 4—2—2 
^1%% + %^/l% + %^/l% 








Pw 


= Im 


2— 2| |2 / 
VlV2\V3\ ( 




2 ~\V2 


2 )] + eye. 






Pn 


= Tie 


/ 4—2—2 
{V1V2V3 ~ 


lA^sl 2 ) 


(M 2 


+ \V3 


I 2 ) 


+ eye. 


P12 


= Im 


4—2—2 (\ |2 1 1 |2\ 

V1V2V3 (ml + M J 


+ eye. 






P13 


= Im 


2— 2| |4 / 


J»7i| 


2 ~\V2 


2 )] + 


- eye. 







(5.15) 
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where "eye." denotes additional terms with the indices cyclically permuted and we have 
noted that four of the generators have extra continuous symmetries. Also, P5, P9, P\\ and 
P12 have a Z 6 symmetry in addition to the requisite U(l) x Oh for fixed r = P^ 2 . An 
arbitrary Oh(Tj 5 ) invariant polynomial may be expressed in terms of the generators (|5.15|) 

P(Vi,V2,Vs)= E C^l n5 Pifl n5 ( Vl , m , m ) (5.16) 

ni,...,n 5 ,X 

where the basis is given by 

p(X) pni pn 2 pn 3 pn 4 pn 5 p /c 17\ 

r n±,...,n 5 — r l r 2 r Z r 4 r h r X \°- Li ) 

where P x = 1, P 6 , P 6 2 , p ePy, P 6 3 , P7, P%, ■ ■ ■ , Pls- This basis has not been constructed previ- 
ously. 

The free energy for 0/i(Tj 5 ) is expressed in terms of the invariant polynomials (|5.15|) : 



p — V F {x) P {x 

Z-^i ni,...,ri5 m, 



•,n5 
ni,...,ri5,X 

aPi + + (3 2 P 4 + (3 3 P 2 + 7 iP! 3 + 72P1P4 + 73P1P2 + 74P3 + 75^6 + 76^7 + 



(5.U 



where the sum runs over the indices described above. These coefficients agree with those 
used by Sigrist and Ueda J| to fourth order, and they do not consider the free energy for 
rj 5 of Oh at higher order. 

C. Minimization of the Free Energy 

The physical gap function minimizes the free energy. The minimum determines both the 
magnitude and the direction of A in representation space. The magnitude depends on the 
parameters a, (3, 7, ... in a complicated fashion. Fortunately, its exact value is unimportant. 
It is zero above T c , small just below the second order critical point and possibly large at low 
temperatures. The direction in 77 space is more interesting, since it determines the symmetry 
of the superconducting phase. 

Theorem § in Appendix [A] guarantees that regardless of the value of the parameters 
a, /3[, ... at least one pair of critical points of the free energy lies on each rotational symmetry 
axis of the representation space (r] space). The theorem does not say which of these critical 
points if any is the absolute minimum, but sufficiently close to the critical point one of 
them is |p9|| . This is a consequence of the Morse Theory of critical points combined with an 



accounting of the critical points of a fourth order polynomial (P) in terms of point group 
orbits. Even as the magnitude of A grows, an intermediate symmetry phase is the ground 
state for most values of the parameters. 

1. ID Irreps 

Because of its simple form, the minimization of the free energy for one dimensional 
representations ( |5.4j ) is straightforward in principle. To fourth order, the free energy is given 
by 
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F = a\r h \ 2 + p\r ]l \* + --- (5.19) 
For a ~ (T — T c ) < and /3 > 0, the minimum is rjf = r e 2ie with 




r o = J!ri + --- ( 5 - 2 °) 



where 9 is an arbitrary phase angle which parameterizes the ground state degeneracy. This 
breaks the U(l) symmetry to Z2, but it does not break the point group symmetry. 

Higher terms in the free energy may be considered as well. The Ginzburg-Landau for- 
mulation is best near the second order critical point where it is well-known that the fourth 
order free energy (perturbed by a few higher order terms to break the residual degeneracy) 
provides a good description of the system. At lower temperatures, the magnitude of the gap 
grows and higher order terms become important. Eventually the perturbative expansions 
in A/uj c and Vzz, may break down due to a finite radius of convergence and asymptoticity, 
respectively. Also, the exact form of the temperature dependence of the coefficients becomes 
important, so it no longer suffices to make the ansatz that the coefficients are independent of 
temperature except for a ~ (T — T c ). Nevertheless, there does exist an effective free energy 
even at low temperatures which is related to the perturbative Ginzburg-Landau free energy 
through resummation, and information about the system at low temperatures (especially 
the symmetries) can be extracted from the higher order terms. 

As these terms are considered, the magnitude of the gap function ro(a, /3, . . .) takes values 
on a branched cover of the parameter space. Consider the minimization of the sixth order 
free energy, F = a|^i| 2 + /3|?]i| 4 + 7|?7i| 6 + -- -, which is of interest in the case of higher 
dimensional irreps. The minimum becomes 



r = ^ (-sgn(/3) + V/1 + 3M7//3 2 ) ^ + • • ■ (5.21) 

The value of the free energy at the minimum is conveniently expressed in terms of the 
function 

*(x,y) = — P^{2(l + 3y) 3 / 2 -sgn0r)(2 + 9y)} 

2 2[ X \V L J (5.22) 

= -^(i-ly + TeV 2 ----) x>0,\y\<l 

where the series is its critical (small a) expansion. The free energy is given by 

F mm = a 2 $([3,\a\ 1 /{3 2 ). (5.23) 

The function $>(x,y) is convenient because it is a monotonically increasing function of both 
\x\ and y. 

It is possible to have first order phase transitions in this parameter space at eighth order, 
where a small change in one parameter causes a large change in the gap function because 
the minimum hops from one sheet to another. The location of these critical points is non- 
universal. The symmetry of the superconducting state need not change at this type of phase 
transition, and we do not consider them further. Note that away from the transitions, the 
higher order terms simply renormalize the leading non-zero coupling f3 according to ( 5.22|) , 
and we define f3 R = -l/(4<J>(/3, \a\j/(3 2 )) so that F min = -a 2 /(A(3 R ). 
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2. 2D Irreps 



The free energy for the two dimensional irreps has the same general features, but it also 
allows the point group symmetry to be broken. We express F in terms of one real and one 
complex variable 

r 2 = \V[\ 2 + W 2 \ 2 = % (5-24) 

and conversely, 

v [ = re 2id+w/2 sech l/2 (Tie w) /y/2 r( 2 = re 2ie - w/2 sech 1/2 {Tie w) fy/2. (5.25) 

The remaining degree of freedom in rj[ and r)' 2 is a phase, 6. This is the zero mode of the 
Goldstone boson which is irrelevant by gauge invariance. These variables behave nicely under 
Oh x U(l). For example, in of Oh, C^ a : w — > w + 2m/ 3 and C^ z : w — > —w. Consider 
the two dimensional irreps rj of Oh and T^^ of D&h first. The generators of the invariant 
polynomials for these irreps (|5.8|) become 



P l = r 2 

Pi = r 4 sech 2 (IZe w) 

(5.26) 

P 3 = r 6 cos(Tm3w) sech. 3 (lZew) 

P4 = r 8 sm(Tm 3w) sech 3 (IZe w) taiah(TZe w) 

which are manifestly invariant under the symmetry operations. In the (r, w) variables, the 
free energy is 

F = ar 2 + fir* + (3' 2 r%ech 2 (Tie w) 

12 v ' (5.27) 

+7i r 6 + 72 r 6 sech 2 (7?.e w) + 73 r 6 cos(Xm 3w) sech 3 (7?.e w) + ■ ■ ■ 

with the same coefficients defined in Eq. (|5.1C| ). 

The order parameter A(r, w) is non-zero below the critical point, breaking the U(l) sym- 
metry. The value of w (the direction of A in rj space) below T c determines whether and how 
the point group is broken. A glance at the functions P\, . . .P4 ( |5.26| ) which generate the 
invariant polynomials reveals that the generic extrema are w = +00, —00, 0, in/3, . . . , bm/3. 
These are the critical points identified by Theorem 2 in Appendix [A], since the rotational sym- 
metry axes in the two dimensional rj representation space are (1, 0), (0, 1), (1, ±1), (1, ±u) 
and (1, ±u 2 ). They give a non-trivial residual symmetry. 

The minimization of the free energy is straightforward. We consider F to sixth order. 
The imaginary part of w only enters through cos(Xm3w), so 

^Httf z r «>o (5 - 28) 
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where the renormalized coupling 73^ = 73 at sixth order. At this order of perturbation 
theory, Xmw is unconstrained if 73 = 0. The free energy takes the form F = ar 2 + 
(3[r A + /3 2 r 4 :r 2 + 71 r 6 + 7 2 r 6 x 2 — |7 3 | r 6 x 3 to sixth order with Xmw given by ( |5.28| ) and 
x = sech(TZew) G [0, 1]. The minima either have x = or x = 1, since the coefficient of 
x 3 is negative and there is no term proportional to x. They correspond to phases in which 
the point group is only partially broken, and A points along a rotational symmetry axis in 
77 space. Specifically, these phases are 



x = 



A = r k x k y k z (k z + ujk x + uj k ) 

A = r k x kyk z (k 2 z +u 2 k 2 x + ujk 2 ) tOT l3 ° tUh 

A = r (k x ±ik y ) for of D 6h 



A = r k y k z (k x ±ik y )(k 2 -3k 2 x ) for Tg of D 



6h 



A = r k x k y k z (2k 2 z - k 2 x - k 2 y ), eye. (73 < 0) 
A = V3r k x k y k z (k 2 x - k 2 y ),cyc. (73 > 0) 

A = r k Xl r {\k x ±^ky) (73 <0) 



(5.29) 



T , 7^ , \ / \ iur 1 t ui i^6h 



A = rofey.ro^i^) (73 > 0) 

\> x kyk z (ky 3k 
A = r Q k 2 k z {k 2 -3k 2 x ),--- (73 >0) 



for T 3 of Ot 



for I\ of Dp 



A = nW^(*J-3«),- (73 <0) forr - ofD6h 



where r is the value of r for the ground state, and u cyc." means terms of the same form 
with the indices cyclically permuted. 

The magnitude of the ground state gap function, r , is given by an expression of the same 
form as Eq. ( |5.21| ), with /3 = f3[ and 7 = 71 for x = and f3 = /3[ + (3' 2 and 7 = 71 +72 — I73I, 



for x — 1. It is natural to define the following renormalized couplings 



PlR 4$(/?{,M7iM 2 ) 

P2R = —ftlR + 



-1 



(5.30) 



4$(# + ^, |a|( 7 i + 72 - |73fl|)/(/3i + /? 2 j 



»2^ 



so that the free energy of the ground state is given by F min {x = 0) = — a 2 /(4/^) and 
F m in(x — 1) — — a V(4(/?i + /3 2 ))- The ground state is x = when (3 2 r > and x = 1 when 



/?2R < 0. The ground state gap functions for the three distinct phases are shown in Table [If I 
Typically, the value of r changes discontinuously at the boundary between these phases, so 
the transition is first order, as expected when the symmetry does not break to a subgroup. 

At eighth and higher orders, there are small regions of the parameter space in which 
the ground state has < x < 1 and the point group is broken completely: T^(Oh) — > 
Ti (D 2 h),^5 (D 6h ) — > r^(Cj) and ^(Dq^) — > Ff (C*), where Cj is inversion. Second order 
transitions to these phases can result from frustration due to the competition between the 
terms minimized at x = and those minimized at x — 1 (or from the contribution of the 
symmetry breaking generator P 4 ). Consider the eighth order free energy. It has the terms 
f3' 2 r 4 x 2 and ^r 8 x 4 which compete when j3 2 < and £4 > 0. The free energy is a fourth degree 
polynomial in x, which decreases as x increases from 0, then reaching a minimum, it increases. 
If the minimum occurs at x > 1, then the ground state is x = 1 since x = sech7?.ew < 1. 
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On the other hand, if the minimum occurs at x < 1, the ground state is the low symmetry 
phase. We can estimate the critical point taking \(3 2 r\ -C faR, for which the low symmetry 
phase exists when 

W >(«h)" 2 (5.31) 

with fa R < and <5 4 > 0. Since a ~ (T — T c )/T c is small near the initial superconducting 
critical point, this secondary phase transition would occur at much lower temperatures. 

The other two dimensional irrep is Tg of D±h- Again, the simple groups act nicely on the 
projective variables, C 2y : w — > w + in and C 2e '■ w — > —w. The rotational symmetry axes 
in the two dimensional representation space are (1, 0), (0, 1), (1, ±1) and (1, ±z), which cor- 
respond to w = oo, — oo, 0, in, in/2 and —in/2. The generators of the invariant polynomials 
are 

P, = r 2 

P 2 = — r 4 sin 2 (Xm w) sech. 2 (1Zew) 

y ' K 1 (5.32) 

P 3 = r A sech 2 (IZe w) 

P4 = r 6 sm(Im 2w) sech 2 (7Ze w) tanh(7?.ew) 
which are manifestly invariant. The free energy becomes 

F = a r 2 + (3\ r 4 — (3 2 r 4 sin 2 (Imw) sech 2 (Tie w) + far 4 sech 2 (IZew) 

+71 r 6 — 72 r 6 sin 2 (lmw) sech 2 (Hew) + 73 r 6 sech 2 (Hew) (5.33) 
+74 r 6 t&iih(7lew) sech 2 (Hew) sm(Im2w) + ■ ■ ■ 

The structure of the solution is more complicated than in the case of the other two di- 
mensional representations because the generator P 4 that explicitly breaks the intermediate 
symmetries occurs at sixth order rather than eighth order. This leads to low temperature 
phases in which the point group is broken to I. 

We restrict our attention to the sixth order free energy with no symmetry breaking 
parameter: 74 = 0. The relative phase for the ground state is given by 

Imw = {J1 f R <l (5.34) 

with [3 2 r = fa + ■ ■ • (See Eq. ( [5.35|) ). The free energy reduces to FiXmw = n/2) = 
ar 2 + /?i r 4 — (3 2 r A x 2 + fa r 4 x 2 + 71 r 6 — 72 r 6 x 2 — 73 r 6 x 2 , and FiXmw = 0) given by the 
same expression with fa = 72 = 0. Note that the free energy is quadratic in x. The phase 
boundaries are again best expressed in terms of renormalized couplings, 

^ = 4$(A, H71/A 2 ) 

fan = -fan + Am + Mal{fyi + fy3)m + m (5-35) 



faR = ftlR + faR + 



4$(A -fa + fa, |a|( 7l - 72 + 7a)/(A ~fa + fa) 2 ) 
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to sixth order, so that F min = — a 2 / (4/3) with j3 = j3i R for w = ±00, (3 = (3\ R + (3 3R for 
w = 0, *7r, (3 = Pir — P2R + P3R for w = ±in/2. The ground states are w = 0,iir when 
@2R, Pm < 0, w = ±00 when /3 3R > max(0, fan), and w = ±in/2 when (3 2 r > max(0, (3 3R ). 
Expressions for the gap function in terms of the original variables for each of the three phases 
are given in Table fTJ. 



3. 3D Irreps 

In analogy with the analysis of the two dimension irreps, the free energy for the three 
dimensional irreps and T§ of Oh should be expressed in terms of complex projective 
coordinates. There are two relative magnitudes and two relative phases, so 771,772 and 773 
should be expressed in terms of a two dimensional vector in complex projective space w G 
CP 2 ; however, there is no known vector that behaves simply under the Oh operations. The 
best we can do is to use the homogeneous coordinates, i]j/(re 2ie ), and there is no apparent 
simplification of the free energy. 

The ground state of the free energy to sixth order is given by a generic critical point 
provided the symmetry breaking parameter 76 is taken to be zero. The generic critical points 
are (0, 0, 1), (1, 1, 0), (1, 1, 1), (1, i, 0), (1,uj,u> 2 ) and (l,u,u 2 ) up to symmetry, according to 
Theorem 2 of Appendix [A[ At these points, the invariant polynomial generators ( |5.15| ) have 
the following values 



77 


Pi 


Pi 


P 3 


Pi 


Pb 


p 6 


-P7-13 


(0,0,1) 


1 








1 











(1,1,0) 


1 


1 

4 





1 











(1,1,1) 


1 


1 

3 


1 

27 


1 











(l,w,w 2 ) 


1 


1 

3 


1 

27 








1 

6 





(1,7,0) 


1 


1 

4 


















at r = 1. Let (3\ R = — l[4$(/3i, |a|7i//3 2 )] 1 , and define the other renormalized couplings 
such that the free energy of each phase is given by 

^(0,0,1) 



-a 2 








-a 2 


4(/9ik 


+ 




-a 2 




+ 02R + 




-a 2 



^KHiR ~r i^Rj 

= Jm Zl* . i« S + ( o7« Z~a — 1 1/3 S ) ( 5 - 36 ) 



c, -« 2 fliR 757? \ 



a 



4(/5ia + \fon) V 27 6 J \2{p 1R + l(3 3R ) 
-a 2 



F(M ' 0) 4(P m + P2R + |Am) + " ' 

where the couplings are defined sequentially, /3 2 r by F(ofi,i), then j3 3R by i<(i,i,o), then 7 4R 
by -^(1,1,1) an< l finally 75/?, by Fn W(tt ,2). The free energies may be expressed in terms of the 
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original (bare) couplings through the $ function as before, using the coefficients listed above. 
Note that (3jR = (3j to fourth order, but the renormalized 7's receive corrections even at sixth 

order, j m = 74 H and 75^ = 75 H . 

The ground state of the cubic system is shown in Table |II|, and the phase diagram is 
shown in Fig. f| taking 74^ = 75^ = 7 6 # = 0. When 74^ and 75^ are non-zero, the phase 
boundaries shift by a small amount. They also have a more interesting qualitative effect in 
that they allow the possibility of stabilizing (1, 1, 0) as the ground state. At fourth order, even 
though (1, 1, 0) is a rotational symmetry axis, the (1, 1, 0) phase only exists on the negative 
(3 2 r axis (see Fig. |j) where it is degenerate with the (0, 0, 1) and the (1, 1, 1) phases. This 
degeneracy is the result of the enhanced 5*0(3) x U{1) symmetry (see Appendix [Cp. The 
fact that the fj = (l,l,0)/v2 is only a saddle point at fourth order and not a minimum 
follows from Morse Theory constraints. This is no longer true in the larger parameter space 
of the sixth order free energy, and the ground state is allowed to be (1, 1, 0). Taking /3 3 = 0, 
we find this occurs for 74 > — ff 73 > 0. 



D. Weak Coupling 

This completes the construction of the generalized phase diagrams for SSS with cubic, 
hexagonal and tetragonal symmetry in the normal phase. We are now in a position to reex- 
amine the superconducting phase transition in the weakly coupled theory. Using standard 
techniques j|0] , we find that the weak coupling expansion of the free energy for A in an irrep 
of G is given by 

1 t 00 

F = ^ + ^(0)log-(|A(^)| 2 )^+5:^ m+2 (|A(,; fc )| 2m+2 )^ 

A I c m -l 



_ 1 (2m -1)11 ( 2m+1 A ci2m + l J-PlY Ni{) , (5 ' 37) 
F 2m+ 2--- j- n , (2 -l]C(2m + l) — j N(0) 



2 (m + 1)! V J *v > ^ 8tt 2 

where Fn is the free energy of the normal phase and (3 C = \jksT c . Note that weak coupling 
means that the parameter 0^.\^ RMS is small. As in Sec. |T|.B, the simplest ansatz is to 
consider an interaction of the form ( |3.12j ) 



V F /6 

where we have rewritten it in terms of the Fermi velocity appropriate for non-spherical Fermi 
surfaces. The coefficients in the free energy ( |5.18|) are then integrals over Fermi surface 
harmonics [^TJ, which are found to equal a = |iV(0) log(T/T c ), /3i = j^F^, j3 2 = 75-^4,71 



^-^6)72 = and /3 3 = 73 = 74 = 75 = assuming a spherical Fermi surface. In this 
case, the free energy reduces to the form invariant under SO(3) x £7(1) which is studied 
in Appendix 0, and the weakly coupled system sits at the point on the positive f3 2 R axis 
near the origin, denoted by an "x" on the phase diagram in Fig. f|. As the Fermi surface is 
deformed outward at the diagonals inducing a positive hexadecapole moment, /3 3 becomes 
positive, and (1, z, 0) becomes the ground state. On the other hand, if the Fermi surface is 
deformed inward, (3% is negative, and (l,u,u 2 ) is the ground state. Therefore, we identify 
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(1, i, 0) and (1, u, u 2 ) in T 4 of Oh as the leading candidates for the SSS ground state in cubic 
systems. In hexagonal and tetragonal systems, the candidates are A oc (1, i, 0) • V}~ in of 
D 6h and D Ah . 

VI. THERMODYNAMIC PROPERTIES 

The gap function may be zero for certain values of k. This can occur accidentally, 
but in some cases the gap function is required to vanish by symmetry. Such zeros are robust 
and have a marked impact on the low temperature behavior of thermodynamic quantities 
such as the heat capacity and acoustic attenuation. The characteristic e~ 2A ^ kT exponential 
behavior due to the finite gap, changes to a power law behavior when the gap function has 
zeros. 

The gap function for SSS always has zeros. This is a consequence of its odd parity 
and the fact that it has a single spin component. The zeros are guaranteed by topological 
considerations very reminiscent of the so-called "Hairy Ball" Theorem (the Hopf or Poincare- 
Hopf Theorem) which relates the total index of a vector field to the Euler characteristic of 
the underlying (closed, orientable) surface [41|. Consider the field e*^ = A/|A| where is 



real, which is a well-defined field on the Fermi surface assuming that A has no zeros. It 
takes values on the unit circle in the complex plane. Since A changes sign under inversion, 
4>{—k) = <p(k) + (2k + l)7r. Following the value of 4>(k) as k is taken around any equatorial 
circle on the (closed) Fermi surface, we find that <fi comes back to itself up to an odd winding 
number, 5(f) = 2ir(2k + 1). The non-zero winding number prevents the equatorial circle from 
being contracted smoothly to a point on the Fermi surface. It must encounter a zero of 
A. Although the gap function may have zeros in singlet and triplet superconductivity, only 
in SSS is it guaranteed to have them. The gapped states of triplet superconductivity (the 



Balian-Werthamer ground state j42|) do not occur because A only has one spin component 
in SSS. 

These generic zeros of the gap function may be found using Theorem [I] in Appendix [X[ 
They are given by fixed points of elements of the residual symmetry group H with a non- 
trivial character, and are tabulated in Table [V]. Isolated point nodes arise as fixed points 
of rotations, whereas lines of zeros are associated with reflections. The nodal structure of 
A determines the density of states near the Fermi surface and consequently the scaling of 
thermodynamic quantities. With a few assumptions the scaling exponents may be computed. 
For example, the heat capacity, which for a constant gap vanishes exponentially, scales as 
T 3 , T 2 and T with point nodes, line nodes and vanishing order parameter, respectively, in 



a defect-free superconductor |27|]. Multiple line nodes lead to TlogT This power law 



scaling is a hallmark of unconventional superconductivity |27j. It must occur in SSS. For 
example, the candidate ground states (1, i, 0) and (1, u, u 2 ) of of Oh discussed above have 
point nodes, so their intrinsic heat capacities scale as T 3 . 

VII. SUMMARY 

The occurrence of a HM AFM normal state, with one conducting and one insulating 
spin channel, has been shown to provide the possibility of a novel superconducting phase for 
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which the operation of time reversal has no part. The state is best considered as a condensed 
phase of spinless fermions, with a gap function that is odd in k and a pair wave function 
that is odd upon interchange of particles. The form of anisotropic interaction required to 
form this state was obtained, and the resulting gap equation is found to be of the BCS form. 
The allowed symmetries of the gap function have been enumerated for cubic, hexagonal and 
tetragonal lattices, and the corresponding conditions on the parameters of the free energy 
have been determined. The quasiparticle spectrum is necessarily gapless. For point nodes or 
a line of nodes, this gaplessness gives rise to power-law behavior in T or uj. For intersecting 
lines of nodes, there will be logarithmic terms such as T log T as enumerated by Nazarenko 



It is anticipated that examples of such phases can be found in transition metal com- 
pounds. The recent suggestions that Sr 2 Ru04 may be displaying triplet superconductivity 
|44| , together with predictions of HM AFM states in transition metal oxides ||19|| , both in- 



dicate that transition metal oxide compounds present a favorable possibility of obtaining 
single spin superconductors. 
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APPENDIX A: TWO THEOREMS 



The minimum of the free energy determines the order parameter of the superconducting 
phase, and the zeros of the gap function in turn determine the scaling of thermodynamic 
properties with temperature. The generic values of these minima and zeros are fixed by 
the symmetries of the system, and they may be determined without resorting to explicit 
representatives of the symmetry. 

This Appendix presents two theorems which are useful in this regard. Theorem [I] may be 
used to find the generic zeros of the gap function. The zeros arise as fixed points of elements 
of the residual symmetry group H that have non-trivial character. This is a refinement of 
the procedure used by Volovik and Gor'kov ||, who identified the zeros with specific group 
elements. The utility of our theorem is that the appropriate group elements may be read off 
of standard character tables. 

Theorem ^| may be used to find the direction of the generic critical points of the free energy 
in representation space. Note that the magnitude of the solution (the magnitude of the gap 
function) is not determined, but this does not affect the symmetry of the superconducting 
phase. The theorem also applies to more traditional applications, such as to show that the 
Fermi surface is orthogonal to the axes of rotation of the crystal where they intersect. This 
theorem applies to the full non-perturbative free energy. 

Theorem 1 Suppose A (A;) is in the irrep V of the little group H , and suppose g G H has a 
non-trivial character, Xr>{g) ^ 1- Then A(ko) = for any fixed point ko of g (i.e. for any 
ko such that g ko = ko). 
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Proof: Since H is the little group, diml" = 1 and 



g : A(fc„) - A(gk Q ) = X r(g) A(fc„) (A.l) 
But k is a fixed point, so 

(i-xAg))A(k ) = o. (A.2) 

And we arrive at the result, A(k ) = provided Xr'(g) 7^ 1- 1=1 

This theorem provides a relatively easy means to identify the nodes of the gap function 
guaranteed by symmetry. As an example, consider the (l,u,u 2 ) state of of Oh, which 
transforms as of the little group (see Table The characters of the elements of C^i 
are listed in common character tables |30|. C3, C^" 1 , /, Sq and S*^" 1 have non-trivial characters 



and of those, C3 and C^ 1 have fixed points, the two points where the axis of rotation (1,1,1) 
intersects the Fermi surface. 

This information is tabulated in Table [V]. Note that in some cases any function in the 
specified irrep must have a divisor whose little group is larger than the little group of the 
function itself. Theorem [I] applies to these factors as well, and any zero of the factor is also a 
zero of the function. An example is kyk z (ky — 3/c 2 .) which transforms as Tg of D 6 h- Its little 
group is D 2 h, but it has the factor k y k z (ky — 3k^) that transforms as the one dimensional 
irrep of D^. This accounts for two additional lines of nodes. 

Theorem 2 Suppose that rj transforms as the irrep F ofG, a subgroup of Oh or D^h, where 
7] is a dimT dimensional complex vector and suppose that the U(l) invariant function F(rj) 
is in the trivial irrep Tf of G. Also suppose that there exists an element g G G with a fixed 
point r] up to a phase: gr]o = e^rjQ. Then e^g lies in the little group of the gradient . 

Note that dimT' < 3. Also note that J^- = djF is the complex gradient, which could 
be expressed in terms of its 2 dim T real components. The theorem states that all of these 
components are invariant under e^g where the phase is equal to the one that appears in the 
i] transformation. 

Proof: 

9 : djFivo) -flfFfo)^ = g^djFivo) (A.3) 
Since gr] = e^rjo, we have 

e i *g jk d k F{r jo )=d j F{ no ) (A.4) 

Thus, djF(rjo) is invariant under e % ^g\ that is, e l ^g is in the little group of djF(rj ). □ 

This theorem is particularly useful when g is a non-trivial rotation (and 7^ 7r); i.e. when 
det = +1 and g is not the identity matrix. Then the fixed points of g are on the axis of 
rotation. In this case, the theorem states that on an axis of rotation, djF points along the 
axis: djF(r] ) oc fj . If we restrict to a surface of constant \r]\ 2 = r 2 , then we find 



OF 



= (A.5) 

\v\=r 
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on an axis of rotation. The axes of rotation are critical points of F in the angular variables. 

We again consider (1, u, uj 2 ) of of Oh as an example. This is a fixed point of C 3 about 
the (1,1,1) axis up to a phase which is a third root of unity: C 3 (l,u,u 2 ) = u(l,u,u 2 ). 
According to the theorem, the gradient of the free energy at (1,u,uj 2 ) must be invariant 
under uC 3 ; hence, it is proportional to the dual vector (1, u 2 , u). This radial vector projects 
to when r is held fixed, so (1,uj,uj 2 ) is a critical point of the free energy. 

Note that perturbation theory is not used directly in the proof of these theorems. To the 
extent that we can say that the functions of interest lie in a particular irrep, the theorems 
are non-perturbative. 

Theorem |2| also has implications for solutions of the gap equation. Since it is an auxiliary 
equation of motion for the free energy, the theorem implies that the full gap equation is 
stable for A (or more precisely, r] m (r)) pointing along any of the rotational symmetry axes 
of the representation space. Stability means that the function 

-^E^^tanh(i^) (A.6) 

is invariant under Q. This reduces the gap equation to a one dimensional problem, only 
slightly less tractable that the usual singlet BCS solution. 



APPENDIX B: SOME EXACT RESULTS 

In this Appendix we present an exact solution of the gap equation at zero temperature for 



two of the cases studied in Section |T|. These two gap functions, d = (0, 0, 1) and d = (1, i, 0) 
of r^" of Oh, are plotted in Figure |3], along with the IY gap function, which has not been 
calculated in closed form. 



1 _ n „„n„\,(„U/o\ , l„./q , ,,9 I V 



3 



arcsinh(z//v / 3) + £zV3 + v 2 + ;r7 = log ( —f=- — /n 2 J (B.l 



A ( o,o,i) 6 6^ Vv / 3 + v / 3 + ^ 

where v = u c /A rms . 

1 v 2 v(9-2v 2 ) (13 



It is evident from the plot that the magnitude of the gap function at low temperature 
decreases as the number of gap function nodes increases, but that the effect is not dramatic. 



APPENDIX C: THE ROTATIONALLY INVARIANT FREE ENERGY 

In the text we constructed and analyzed the free energy with the gap function in various 
irreducible representations of the cubic, hexagonal and tetragonal point groups. The result- 
ing phase diagram is quite complicated because of the many couplings necessary to specify 
the pairing interaction. Only a small subset of these couplings corresponds to physical per- 
turbations. 
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The free energy simplifies greatly if we make the natural ansatz that the dominant channel 
for pairing comes from the interaction ( 3.12 ) 



k ■ k' 



w^^-Mt^, (c.i; 



and that the Fermi surface is spherical. This interaction has an 0(3) symmetry, a higher 
symmetry than we have considered in Section |V[ Let the gap function be in of Oh, the 
only cubic irrep that pairs under the interaction ( |C.1| ). It has two non-cyclic degrees of 
freedom, the magnitude, r, and the angle, ip, between its real and imaginary parts: 



A = r (1 + i cos tp,i sin tp, 0) • k/y/2. (C.2) 

The overall phase and the three Euler angles do not affect the energy. The free energy is 
expressed in terms of the well-known 0(3) invariant polynomials 

r> I 1 2 i I 1 2 i I 1 2 2 

Pi = M + \m\ + ml = r 

r> l 2 , 2 , 2 1 2 4 2 v^' / 

"2 = \Vl + V2 + V3\ = r cos f 

and it is given by 

F = ar 2 + /?ir 4 + (3 2 r 4 cos 2 tp + ^r 6 + <y 2 r 6 cos 2 y? H (C.4) 

According to Theorem |2] (extended to 0(3)), the generic critical points are r (0, 0, 1) • k 
and ro(l, i, 0) • kj v2 up to rotation. Note that these have a line of zeros and 2 point zeros, 
respectively. The free energy has non-generic critical points in small regions of parameter 
space at eighth order in A. 

At sixth order, the ground state gap function is given by 

A = r k z (3 2R < 

A = r {k x + ik y )/V2 (3 2R > (C.5) 

where (3 2R = - [4$(/3i + fa, |a|(7i + 72) /(/^l + faf 

and note that (3 2 r = f3 2 at fourth order. The value of r is given by the expression in Eq. 
([5.21]) with p = (3i and 7 = 71 in the first case, and (3 = (3\ + (3 2 and 7 = 71 + 72 in the 
second. 
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TABLES 



TABLE I. Categorization of normal state symmetries, number of allowed broken symmetry 
broken states, and allowed values of S and L for several fermion liquids (S, L values for the BCS 
case refer to cubic crystal symmetry). "HFS" refers to the picture where the spin is frozen into 
the lattice and therefore is not a separate symmetry of the normal state; it is still unclear if this 
picture gives the best description of the heavy fermion superconductors. "No J" indicates "HFS" 
without inversion (see text). Symmetry group notation is given in the text. 

System Normal Broken Symmetries Pairing Type 



3 He 



CxS xT x Jx U(l) 



oo 



S=0, L=even 
S=l, L=odd 
S=0, L=0,2,4,6 
S=l, L=l,3,5,9 
S,L Coupled 



BCS 



G xS xT xlx U(l) 



Finite 



HFS 



g xT xlx U(l) 



Fewer 



Even or Odd Parity 



No 1 

sss 



g xTx U(l) 

g xix u(i) 



Still Fewer 
Fewest 



Impure States 



L=odd 
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TABLE II. Leading basis functions for the irreps of Oh-, -t>6/i and D^. 

Irrep V Basis t/>(T, m; k) 

Octagonal — Oh 



IT 




1; /c) — k x kyk z {k^ ky)(k 2 k 2 ){k 2 


1 2 


vy 1 2 


1- = ka.k„,k. 


■p— 


n 1 3 


1 • M h h h ( Oh^ i-2 7,2 \ 
— h lx h,yh, z \Ll^ z tv x tVy) 




^C 1 3 


*,rv) — y ohu x i^yH z yiv x ivy) 




V>(r 3 - 


1'; fc) = k x k y k z [k 2 z + wAi 2 + w 2 fc 2j 






2'; fc) = k x kyk z (k 2 z + w 2 ^ + W A 2 ) 




V>(r 4 " 


1; k) = k x 




V>(r 4 " 


2; A) = fcj, 






3; k) = A z 




V<(r 5 - 


1;§=M^"^) 




^(r 5 - 


2; £:) = A y (A; z — k x ) 






3; A) = A^A^ — k y ) 



Hexagonal — D§h 







1; k) = k x k y k z (k x — 3A 2 )(A 2 


r 2 




1; A) = fcj; 


r 3 " 




1; k) = ky — 3k x k y 


r 4 " 




1; k) = k x — 3k x k 2 


r 5 " 


V<(r 5 - 
^(r 5 - 


1; A) = A^ 
2; k) = ky 


r 6 


^(r 6 - 


I5 A?) — k x kyk z {ky 3k x ) 
2; A) = kyk z (ky — 3A 2 ) 



Tetragonal — D^h 



rr 




1; k) = k x kyk z (k x — k 2 


r 2 " 




1; k) = A z 


r 3 " 




1; A) = k x kyk z 


r 4 " 


V>(r 4 


1; A) = (A; 2 , — ky)k z 


r 5 " 


^(r 5 - 


1; k) = k x 




^(r 5 - 


2; fc) = Ay 
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TABLE III. SSS symmetry breaking phases. 



r 


HQ?) 


Hphys 


ft 


# 










Octahedral 


-Oh 


rr 




o h 




1 


k x k y k z (k x — ky)(ky — k^)(k 2 — k x ) 


r 2 


o h (r 2 ) 


o h 


— 


1 


kx ky k z 




JJ4h{L i ) 




/9o n < 'Yq r? > 


3 


.nil- L U fZ~ 2 h2\ 

V <Jh X K,yh Z {K, x tVy) 




U Ah l 1 3 ) 


n 

JJih 


/3or < 0,7qr < 

M Zl\ N "5 / Oil ^ *-* 


3 


U h h fll-Z 1,2 J,2\ 

h, X IS,yh, Z \Lh, z h, x fay) 




I h{^ 2 ) 


r\ 
<Jh 


H Art VJ 


1 


1^ 1^ 1^ ( 7*2 i , ,7^,2 i , ,2 ]~2\ 
K x K y K z \K z + U)K X + LO Ky) 










1 


K x KyK z \K z + Ul K x + ^Ky) 


r 4 " 


c 3l (r 2 ) 


D 3d 


Bin < < /?9R 
Holt ~ J - HArt 


4 


k z + Lok x + (j 2 /c y 






D M 




4 


fe z + to 2 k x + w/cy 






D M 


B<?R* |3qD < 
HArt i Hort - w 


4 


k z -\- k x -\- ky 




-^4/1 l-l 2 ) 




4/?9R < Bin > 

^Hajx -~ Hort-) Hort ^ u 


3 


h 

h, z 




W/U 1 3 ) 


J->Ah 


< < 4/?9R 
w ^ Ho r\. ^ ^H Art 


3 


kgr | iky 




WM 1 4 ) 


u Ah 




3 






n i"r — ~\ 

^2/i (1 2 ) 


L>2h 


/3, R = /?9P < 

Hort VJ 7 H2.it w 


3 


L _l_ 




^2/H 1 4 j 


n , 




3 


L t, 






D M 


/%? < < 


4 








D m 




4 


^z(,^x ^y^ ^ ^xi^ky ^z) ^^y^^z ^x 






Dm 


< 


4 


kz{kx ~ ky) + k x (ky — k^) + ky(k^ — k%) 




£>4ft(rr) 


D 4h 


4/?2R < /3 3 r,/3 3 r > 


3 


kz{k x — ky) 




C4/i( r 4 ) 


D 4h 


< /3 3 R < 4/? 2 R 


3 


k x (k y — k z ) + ik y (k z — k x ) 




cwr^ ) 


D 4h 




3 


kx(k y — k z ) — ik y (k z — k x ) 




£> 2 ft(rj) 


D 2h 


/3 3 i? = 0,/? 2 R<0 


3 


k x (k y — k z ) + k y (k z — k x ) 




Ay,(ri) 


D 2h 




3 


kx(ky — k z ) — ky(k z — k x ) 
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TABLE IV. SSS symmetry breaking phases (cont.). 



Hexagonal — 



IT 




D 6h 




1 


7 7 / ( 1 2 7 2 \ / 7 2 7 2 \ 
K x K y K z {K x — 6K y ){K y — 6K X ) 


17 




D 6h 


— 


1 


k z 


r 3 " 




D 6h 


— 


1 


1,3 r>h,2h. 


1 4 


Dr; (Y~\ 






1 

± 


1,3 qi, 1,2 

"> x 0hj xr"y 




D 2h {T~ 4 ) 




02R < 0, 73/? < 


3 


k x 




D 2h (r 2 ) 


D 2h 


02R < 0, 73R > 


3 


ky 




Cm (Ft) 

^bh\^ 5 ) 


1 -'<ah 


Bon > 


1 


k~. A- ik„ 




CWitTg ) 






1 


k x iky 


r 6 


£>2h(rr) 


D 2h 


&2R < 0, 73R < 


3 


k x kyk z {ky — 3k x ) 




^2h(ri") 


D 2h 


02R < 0, 73R > 


3 


KyK z {Ky — OK x ) 




Cfi/, (To ) 


D 6h 


02R > 


1 


kukrik^r ~\~ zkii}(k„. 3/c~> ) 

y & \ j> y / v y x / 




C6/i( r 3 ) 


D 6h 




1 


kyk z (k x iky){ky 3/c^) 








Tetragonal — 






rr 


£>4fc(rr) 


D 4h 




1 


/ 7 / / / 2 / 2 \ 

K x KyK z {K x — Ky) 




^(r 2 -) 


D 4h 


— 


1 


k z 


r s 




D 4h 





1 


k x ky k z 


r 4 " 


£>4h(rj) 


D 4h 




1 


(k x — ky)k z 




cWrj ) 


D 4h 


02R > 03R, 02R > 


1 


kj£ | %J\iy 




c^rj ) 


D 4h 




1 


kx iky 




D2h(T 2 ) 


D 2h 


02R,PsR < 


1 


kx H~ ky 




D 2h {T- A ) 


D 2h 




1 


kx ky 




D2h(T 2 ) 


D 2h 


PsR > @2R, @3R > 


1 


ky 




D 2h (T A ) 


D 2h 




1 


k x 
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TABLE V. Zeros of the gap function guaranteed by symmetry. 

T H(T') Non-trivial Mappings Generic Zeros 

Octahedral - Oh 

IY O h (T^) I,8S 6 ,3(T h ,S4,6<T d 9 circles: {h = 0} i= i j2 ,3, {h = ±kj} if Lj 

r 2 O h (T 2 ) 6C 4 ,6C 2 ,/,8S' 6 ,3o- fe 3 circles: {/cj = 0}i = i i2 ,3 

D Ah {T^) I,2S4,a h ,2a v ,2a d 5 circles: {h = 0} i= i j2 ,3, {fa = ±k 2 } 

D 4h (Tz) 2C 4 ,2C 2 ",I,a h ,2a v 3 circles and 8 points: {fa = 0}i=i >2 ,3, {(±l,±l,±l)/>/3} 

T h {T 2 ) AC z ACz l J,Za h AS & ASQ l 3 circles and 8 points: {fa = 0} l=1 , 2 , 3 , {(±1, ±1, ±l)/\/3} 

T h {Tl ) 4C 3 ,4C^ 1 ,I,3a h ,4S 6 ,4S^ 1 3 circles and 8 points: {fa = 0} l=1 , 2 , 3 , {(±1, ±1, ±l)/\/3} 

rj C 3l (r 2 ) C 3 , C3- 1 , /, S 6 , S 6 l 2 points: {±(1, 1, 1)/V3} 

C 3 i(r 3 ) C 3 , C3- 1 , /, 5 6 , 5^ 2 points: {±(1, 1, 1)^3} 

D 3d (T 2 ) 3C 2 , /, 2S 6 (<7 ft of A> d ) 1 circle: {fa + k y + fa = 0} 

£>4h(rj) 2C 2 ,2C 2 ',/,2S 4 , ( t /i 1 circle: {& 2 = 0} 

) C4, C 2) C^ 1 , J, 5 4 , 5J 1 2 P° ints: 0. !)} 

C 4h (rj) C 4 , C 2 , Q" 1 , 7, 5 4 , 1 2 points: {±(0, 0, 1)} 

r 5 - C 3i (r 2 ) C 3 , C3- 1 , /, S 6 , S^ 1 2+12 points: {(±1, ±1, ±l)/>/3}, {(0, 0, ±1)} + eye. 

C 3 *(r 3 ) C 3 , C 3 -\ /, 5 6 , 5^ 2+12 points: {(±1, ±1, ±1)/V3}, {(0, 0, ±1)} + eye. 

D 3d (T^) I,2S e ,3a d 3 circles: {fa = fa}^ 

D 4h (T^) 2C 4 ,2C 2 ,I,2S 4 ,a h ,2a d 3 circles: {k z = 0},{fa = ±k y } 

C 4h (Tz) C 4 , C 2 , C 4 \I, 5 4 , 54- 1 2+12 points: {(0, 0, ±1)} + eye, {(±1, ±1, ±l)/>/3} 

C4h(rj) C 4 , C 2 , C4- 1 , /, 5 4 , 54- 1 2+12 points: {(0, 0, ±1)} + eye, {(±1, ±1, ±l)/>/3} 
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TABLE VI. Zeros of the gap function guaranteed by symmetry (cont.). 



Hexagonal — D, 



6h 



IT D 6h (T7 
A^(r 2 

D 2h {T7 

D 2h (F 3 
CW*(r 2 



I, ah, 2S 3 ,2S 6 ,3ad, 3a v 

3C 2 , 3C 2 , 1, ah, 2S 3 , 256 

C 2 , 2Cq, 3C 2 , I, 2Ss,3ad 

C 2 , 2Cq, 3C' 2 ,I, 2Sq, 3a v 

C 2 , C 2 , 1, a' v 
C 2 ,C' 2 ,l,a" u 

r t / t — 1 r t /^—l T C G~l C G~ 1 
^2, ^3; ^3 > 1 > a 3 j <->3 )^6;^6 

I,a v ,a' v ,a"(a' d ,a d of D 6h ) 
C 2 , C 2 , 1, a v (a d , a' d , a" d of D eh ) 
C 3 ,C 3 ,Cq,C 6 ,I,S 3 ,S 3 ,Sq,S 6 ,ah 
C3,C 3 ,C(s,C & ,I,S3,S 3 ,Se,S 6 ,a h 



7 circles: {h = 0}j = i )2)3 , {V3k x = ±k y }, {\/3k y = +k x } 

1 circle: {k z =0} 

3 circles: {k y = 0}, {V3k x = ±k y } 

3 circles: {k x = 0}, {V~3k y = ±k x } 

1 circle: {k y = 0} 

1 circle: {k x = 0} 

2 points: {±(0,0,1)} 
2 points: {±(0,0,1)} 

3±2 circles: {ki = 0}j = i i2)3 + {k y = ±V3k x } 
1+3 circles: {k 3 = 0} + {k 2 = 0}, {k y = +V3k x } 
1+3 circles: {k 3 = 0} + {/c 2 = 0}, {k y = +V3k x } 
1+3 circles: {k 3 = 0} + {/c 2 = 0}, {fcj, = +V3k x } 



Tetragonal — D^h 



IT 


£4/, (IT) 


/, 2S4, a h , 2a v , 2a d 


5 circles: {ki = 0} i= i i2)3 , {k x = +k y } 


r 2 


£>4h(rj) 


2C 2 , 2C 2 , 1, 2S4, ah 


1 circle: {k z = 0} 




^(rj) 


2C4, 2C 2 , 1, ah, 2a v 


3 circles: {ki = 0}j = i j2)3 




£>4h(rj) 


2C4, 2C 2 , /, 2S4,ah, 2a d 


3 circles: {k z = 0},{k x = ±fc y } 




c^OT) 


C4 > C 2 , C 4 1 , /, S4 , S4 1 


2 points: {±(0,0,1)} 




cwrj ) 


C 4 > C2 > C4 , 1, S4, S4 


2 points: {±(0,0,1)} 




£>2 h (rj) 


C 2 , C 2 , 1, a' d 


1 circle: {k x = —k y } 




£> 2h (rj) 


C 2 ,C 2 , 1, a" d 


1 circle: {k x = k y } 




£>2h(rj) 


C 2 , C 2 , 1, a' v 


1 circle: {k y = 0} 




£>2h(ri) 


C 2 ,C' 2 ,I,a" u 


1 circle: {k x = 0} 
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FIGURES 



( a ) Ferro 


: up 




E F Down 




-1.5 -0.5 0.5 1.5 

Energy/Bandwidth 

FIG. 1. (a) Model spectrum with rigid exchange splitting that illustrates a HM FM system, 
(b) Model spectrum for a HM AFM illustrating that the channels must have different structure. 
The peak at the Fermi level is merely an artifact of the form and symmetry of the model. 
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FIG. 2. Total densities of states for each spin direction in the double perovskite compound 
La2VCu06- Top panel: parallel alignment of the Cu 2+ and V 4+ spins, with exchange splitting of 
0.25 eV for the lower lying (Cu) states and 0.5 eV for the higher lying (V) states. Bottom panel: 
antiparallel alignment of the spins, resulting in a HM AFM system with the Fermi level Ep lying 
in the gap of the down spin channel. 
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FIG. 3. The RMS value of the gap function at T = relative to the coupling boson frequency 
cutoff, oo c , for the singlet BCS case, and for three SSS cases. The coupling strength A is defined in 
the text. Note that the point nodes of the "liO" gap, the line nodes of the "001" gap, and the set 
of nine line nodes for the case do not affect the RMS gap value drastically. 
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Phase Diagram 







2R 





r ^ 

F.S. 

v J 




(l,i,0) .... 


(1,1,1) 


(0,0,1) 



3R 



(1,1,0) 



FIG. 4. The Ginzburg-Landau phase diagram for SSS with 5 cubic symmetry has five phases 
near T c . The cubic group is only partly broken in each of these phases. The notation is explained 
fully in the text, but for the gap function is proportional to d- k, where d is the vector labelling 
each phase. The phases at the top of the diagram have point nodes and those at the bottom have 
line nodes. The phase with d = (1, 1,0) is only stable on the negative fan axis (see the text). The 
"x" on the positive fan axis denotes the weak coupling point with a spherical Fermi surface and 



the interaction ( |5.38| ). When the Fermi surface is deformed outward or inward at the diagonals as 
pictured, the ground state is a phase on the right or left side of the diagram, respectively. 
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